The nature of spectral gaps due to pair formation 
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Several phenomenological self-energies have been presented to describe the pseudogap in cuprates. 
Here, we offer a derivation of the self-energy in two dimensions due to pair formation and compare 
it to photoemission data. We then use our results to address several questions of interest, including 
the existence of magneto-oscillations in the presence of the pseudogap, and the two length scale 
nature of vortices in underdoped cuprates. 
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Various models for the self-energy have been presented 
to describe photoemission spectra for the cuprate pseu- 
dogap phase. 1 - The basic functional form is 



E(fc,w) 



oj - X k + iT 



(1) 



where A is the energy gap 2 and Tq the broadening. For 
the pairing scenario, X k = — e^, where e k is the single 
particle dispersion. This Ansatz has a long history go- 
ing back to the original BCS theory^ where it implicitly 
describes broadening due to impurities^ 

Lee, Rice, and Anderson^ were able to derive the same 
functional form for a one-dimensional density wave state, 
with Xf. = £k+Q where Q is the wavevector of the density 
wave. In this one dimensional case, long range order is 
not present (A 2 = (A 2 )). The result was derived at 
lowest order restricting to static thermal fluctuations. In 
this case, T is replaced by T 2 = Vp/£ where Vp is the 
Fermi velocity and £ the correlation length. A similar 
derivation in two dimensions yields instead^ 



ImYi{k, ui) 



A 2 



vV + e^ + r 2 



(2) 



Eq. 1 was proposed some time ago to describe data in 
the pseudogap phase of the cuprates^— In Ref. 0, it was 
motivated by a 'zero dimensional' approximation where 
the fermion dispersion is ignored (i.e., e k -q ~ £fe) when 
doing the momentum integration (E ~ J DG where D is 
the boson propagator and G is the fermion Green's func- 
tion). In this case, Tq in Eq. 1 reduces to that of time 
dependent Ginzburg-Landau theory, and should scale ap- 
proximately as T — T c (as compared to the \JT — T c be- 
havior of r 2 in Eq. 2). This was found to give a good 
account of the T dependence of the photoemission data 
above T c for underdoped cuprates at the antinodal points 
of the Brillouin zone (where the d-wave energy gap is 
largest) X It was claimed in this work that this functional 
form could be motivated in higher dimensions as well, 
but as we show here, this is dependent on the value of 
two physical parameters, A/T c and vp/^qT c . 

Recently, Senthil and LeeiS proposed a related Ansatz 
for the zero temperature limit, which was motivated 
by a desire to address magneto-oscillation data in the 



cuprates. Their Ansatz, though, leads to three spectral 
peaks, as opposed to Eq. 1 that either yields two peaks 
(gapped case) or one peak (gapless case) depending on 
the ratio Tq/A. Their result is similar to a related one de- 
rived for a spin density wave by Kampf and Schrieffer. 11 
We note that both results seem to be at variance with 
the expectation that the energy gap should be confined 
to the ordered and 'renormalized classical' phases, and 
therefore should not be present in the zero temperature 
limit unless ordering is present^ 

In this Rapid Communication, we provide a derivation 
of the fermion self-energy due to pairing including both 
the static thermal fluctuations as in Refs. [5||6| and the 
dynamical fluctuations as in Ref. 17]. Above T c , we find a 
result in two dimensions which contains aspects of both 
Eqs. 1 and 2, with the dynamical broadening (Tq) dom- 
inating over the thermal broadening (T%) if vp/^qA is 
small relative to unity (this ratio is tt in BCS theory), 
where £o is the bare coherence length. With a reasonable 
choice of parameters, we find that it quantitatively fits 
photoemission data for underdoped cuprates. At T = 0, 
we find that for these same parameters, three spectral 
peaks are indeed present in agreement with the work of 
Senthil and Lee, though for BCS parameters, only a sin- 
gle peak occurs. 

To lowest order, the electron self-energy is obtained by 
convolving the pair propagator with the hole propagator: 



E(fc, LO v 



d d 
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( 3 ) 

where D is the pair propagator and Go is the bare Green's 
function (Gq 1 = iuj n — £fc), with the sum over boson 
Matsubara frequencies. In the BCS approximation, D is 
A 2 S(q)S(fl), immediately giving rise to Eq. 1 with To = 
+ . In the absence of long range order, 



D- 1 = N Q (x + $q 2 + a\n m \) 



(4) 



withH x ~ (T - T c )/T c , a - tt/(8T c ) and £ oc v F /T c . 
N is the density of states per unit cell. For T > T c , 
the dominant contribution to the Matsubara sum comes 
from the branch cut of D on the real axis. This leads to 



2 



(coth(fi/2T) - 2T/Q) 

T f d d q/(2Tr) d 



Nq J x + ^q 2 iu} n + i{x + ^q 2 )/a + e q ^ k 



Evaluating (d=2 is assumed from here on), we find 



-iA 2 
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(5) 



(6) 



with A 2 = 2,^2. r 2 = vfVx/Zo and f = 2x/a. Al- 
though this formula does a good job of reproducing the 
filling in of the pseudogap with temperature seen by pho- 
toemission, the T dependence of the spectral gap magni- 
tude is not properly reproduced - in particular, the spec- 
tral peak position exceeds A in magnitude for a large 
range of T. This problem can be traced to the definition 
of A itself. In Ref. 0, this difficulty was avoided in the 
derivation of Eq. 1 by ignoring the q dependence of the 
second term in Eq. 5. By making this approximation, 
this term could be extracted outside the q integral. The 
q integral then reduces to the definition of the fluctua- 
tional gap, (A 2 ). The issue, though, is that the static 
terms giving rise to T2 are ignored in this approximation. 

These troubles ultimately stem from the fact that Eq. 4 
is a low q, low Q approximation of the true pair propaga- 
tor. Use of Eq. 4, though, closely matches the exact result 
if the thermal approximation (coth replaced by 2T/0) is 
used to cut-off the integration, and the q integral is 
cut-off at l/£o-— Evaluating Eq. 5 with the cut-off gives 



T 



In 
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where a = — 1/a 2 , b = —v 2 F /^ + 2i(uj + e k )/a and c = 
(uj + e k ) 2 + xvp/^ 2 . At high frequencies, this reduces to 



■'high 
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(8) 



Noting that in this approximation, the fluctuational gap 
is 



(A 2 ) = -JL^ln- 



1 



(9) 



we now find the proper high frequency behavior of the 
self-energy, (A 2 ) /uj. 

Formally, (A 2 ) has a singular temperature depen- 
dence, but for purposes here, we will simply set its value 
to experiment, noting that photoemission spectra indi- 
cate no temperature dependence of A at the antinode/ 
In Fig. la, the real and imaginary values of the self-energy 
from Eq. 7 versus uj are shown for x = 0.1, and in Fig. lb 
the half width half maximum of the imaginary part is 
shown versus x. The parameters used were A/T c = 4 
and vf/£,qT c = 1. These values were chosen so as to give 
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FIG. 1: (Color online) (a) Self-energy from Eq. 7. Parameters 
are x = (T — T c )/T c =0.1, A/T c = 4 and v F /£ T c = 1. (b) 
Half width of ImE, denoted as F, versus x compared to the 
data of Ref. 0. 
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FIG. 2: (Color online) Spectral functions (e& = 0) using the 
same parameters as Fig. 1 for x — 0.1 and x = 2. 



a good account of the experimental half width versus x 
extracted from fitting photoemission data on underdoped 
cuprates using Eq. l£ We note that the A/T c ratio of 4 
(as compared to the BCS value of 1.76) is a typical value 
observed in cuprates. The value of Vp/£oT c of 1 (as com- 
pared to the BCS value of 1.76tt) acts to emphasize the 
dynamic broadening (T ) over the static broadening (r 2 ). 
The ratio of these two values is only 1 /4 (compared to 
the BCS value of tt) and will have further consequences 
as discussed below. As mentioned before, the resulting 
spectral functions have either two peaks or one peak de- 
pending on the magnitude of the half width relative to 
A, with examples shown in Fig. 2. This crossover over 
from gapped to gapless behavior occurs when the ratio of 
the half width to A is about V2, the same as from Eq. 1. 

In the zero temperature limit, the imaginary part of 
the self-energy is given by 
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-q)ImD(Q,q) (10) 



Without cut-offs, this integral is 
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where x in Eq. 4 is now a tuning parameter besides tem- 
perature - magnetic field, etc., i.e., (H — H C 2)/H C 2 - with 
x = corresponding to the quantum critical point where 
long range order appears^ The result is that -ImE al- 
ways grows with frequency, saturating to a constant as 
u) — > oo. The real part of the self-energy can be obtained 
by numerical Kramers-Kronig)A£ and it is found that the 
resulting spectral function is gapless. The reason is that 
formally, the integrals defining (A 2 ) are divergent, so 
cut-offs must be invoked, this time not only in momen- 
tum, but also in frequency as well. We choose to cut-off 
the q integral at l/£o and the integral at l/a. Reeval- 
uating, we find 



ImT, — 



1 



2^N e 

—c(x) tan -1 



Im / dq x [c(x + 1) tan -1 c(x + 1) 

J q x 

c{x)] 



(12) 



where c(y) 1 = yjy + q x + ia(uj - v F q x /£ ), q 1 = 
max(0, (lu — l/a)^o/vp) and (72 = mm(uj£,a/vF , 1) with 
q now expressed in units of l/£o- Similarly, we find that 



(A 2 ) 



1 , 1 {x + l) 2 + l 

+2tan _1 (a; + 1) - 2 tan _1 (a;)] 



x In ■ 



(13) 



(A 2 ) is then used to set the prefactor in Eq. 12. We can 
now evaluate -ImS by doing one numerical integration. 
We note that in this approximation, -ImE vanishes be- 
yond a frequency w c = l/a + vp /£o due to the cut-off in 
il. In fact, we note that the various cut-offs define two 
other frequency scales as well, uj\ = l/a and W2 = vf/£oj 
with lo c being their sum. uj\ is associated with the dy- 
namic part of the pair propagator, and W2 with the static 
part. 

In Fig. 3a, we plot the self-energy from Eq. 12, and 
in Fig. 3b the resulting spectral function, for the same 
parameters as in Fig. la. One clearly see the existence of 
three spectral peaks. We can contrast this behavior with 
that in Fig. 4, where we show the same as Fig. 3, but now 
for BCS parameters. In the latter case, the asymptotics 
of the self-energy sets in at a frequency beyond A, and 
therefore no spectral gap emerges. Similar results are ob- 
tained if one replaces the propagator in Eq. 4 by that in a 
magnetic field in the lowest Landau level approximation. 

Our T=0 results can be compared to the recent work of 
Senthil and Lee^ where a separable approximation for 
the propagator was used. In their work, a propagating 
form was considered 



ImD 



( g 2 + £-3)3/2 



(s(r - n) - 8(t + n)) (14) 



The resulting self-energy at T = is equivalent to that 
for electrons coupled to an Einstein mode with frequency 
r£ That is (w > 0) 
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(15) 
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FIG. 3: (Color online) (a) Self-energy from Eq. 12. Param- 
eters are x = 0.1, A/T c = 4 and vf/£oT c = 1. (b) Spectral 
function (e^ = 0), where a constant 0.1A has been added to 
-ImE so as to resolve the delta functions. 
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FIG. 4: (Color online) Same as Fig. 3, but for A/T c = 1.76 
and vf/^oTc = 1.76-7T. 



where O is the step function. This has a gap between — T 
and +r (with the real part of the self-energy diverging 
logarithmically at ±r). As a consequence, the spectral 
function consists of incoherent peaks at \oj\ > T, and a 
quasiparticle pole within this gap. 

A similar result occurs if one assumes a diffusive be- 
havior which is more appropriate for the disordered phase 



ImD = -- 



2A 2 7r£- 
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The resulting self-energy at T — is (e k = 0) 



(16) 
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where we have used that T = vf/S,- This functional form 
(Fig. 5a) also leads to a spectral function with three peaks 
(Fig. 5b). 

It is interesting to note that in the Senthil and Lee 
formalism, the only energy scale is T, and therefore a 
spectral gap occurs as long as the ratio of T to A is not 
too large. One reason for the difference from our work is 
that in their separable approximation, T is independent 
of q, whereas from Eq. 4, one finds that the relaxational 
rate is strongly q dependent, that is T q — a^ 1 {x+^q 2 )^ 
We also note that formally, the fi integral of Eq. 16 is 
logarithmically divergent when used to define (A 2 ), but 
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FIG. 5: (Color online) (a) Self-energy from Eq. 17, with T = 
0.1 A. (b) Spectral function (ek = 0), where a constant 0.1 A 
has been added to -ImE so as to resolve the quasiparticle pole. 

when calculating the self-energy, this is compensated for 
by the convergence of the q integral in this separable 
approximation. That is, Eq. 17 is well behaved without 
the need to explicitly invoke cut-offs. 

We now turn to the question of the electron pockets 
observed by quantum oscillation experiments.— How can 
such pockets survive in the presence of a large pseudo- 
gap, since these electron pockets should originate in the 
antinodal regions of the zone?— As Senthil and Lee point 
outfit as one indeed finds a central peak inside the gap 
in the low temperature limit, the existence of magneto- 
oscillations is not a surprise (though in our case, we find 
a spectral gap only if the asymptotics of the self-energy 
sets in below A). More generally, quantum oscillations 
are seen in type II superconductors, sometimes for fields 
much less than H C 2- At a semiclassical level, this can be 
understood since the expectation value of the supercon- 
ducting order parameter averages to zero over a cyclotron 



orbit due to phase winding around the vortices. As a 
consequence, type II superconductors are gapless at high 
magnetic fields, with the energy gap causing a broadening 
of the Landau levels. Quantum mechanical simulations 
have demonstrated the evolution of the low energy vor- 
tex core bound states into Landau levels as the field is 
increased^ and similar calculations have been used to 
address the quantum oscillation data in the cuprates. 21 
Extension of these methodologies to a potential vortex 
liquid phase above the resistive H C 2 would be illuminat- 
ing. We remark that the gapless peak in our work (and 
Senthil and Lee's) traces out a large Fermi surface, and 
therefore density wave formation would have to be in- 
voked to explain the small electron pockets that are ac- 
tually observed^ 

We note that for the parameters in Figs. 1-3, the value 
of vp/irA is smaller than £ by a factor of 4ir. If we iden- 
tify the former with the size of the vortex core and assume 
a typical value of 30 A, then the latter is approximately 
400 A. Such a long length has been identified from ter- 
ahertz conductivity measurements j22 and implies a large 
'halo' which exists around the vortex cores, leading to 
the concept of cheap, fast vortices, with the resistive H C 2 
where these halos overlap i 10 ' 24 Therefore, a large A/T c 
ratio and a small vf/^qA ratio are conducive to obtain 
an extended regime above T c and H C 2 where an energy 
gap exists without long range order, a regime that should 
be characterized by fluctuating vortices. 

Work was supported by the U.S. DOE, Office of Sci- 
ence, under Contract No. DE-AC02-06CH11357. We 
thank Mohit Randeria, Todadri Senthil and Patrick Lee 
for discussions. 



1 M. R. Norman, A. Kanigel, M. Randeria, U. Chatterjee 
and J. C. Campuzano, Phys. Rev. B 76, 174501 (2007). 

2 A implicitly depends on k, i.e. A*,. This is allowed under 
the assumption that the pair wavefunction factorizes as a 
product of relative and center of mass degrees of freedom. 

3 J. R. Schrieffer, Theory of Superconductivity (W. A. Ben- 
jamin, New York, 1964). 

4 A. A. Abrikosov and L. P. Gor'kov, J. Exp. Theor. Phys. 
12, 1243 (1961). 

5 P. A. Lee, T. M. Rice and P. W. Anderson, Phys. Rev. 
Lett. 31, 462 (1973). 

6 K. Maki, Physica C 282-287, 1839 (1997). 

7 M. R. Norman, M. Randeria, H. Ding and J. C. Cam- 
puzano, Phys. Rev. B 57, R11093 (1998). 

8 M. Franz and A. J. Millis, Phys. Rev. B 58, 14572 (1998). 

9 J. Maly, B. Janko and K. Levin, Phys. Rev. B 59, 1354 
(1999). 

10 T. Senthil and P. A. Lee, Phys. Rev. B 79, 245116 (2009). 

11 A. P. Kampf and J. R. Schrieffer, Phys. Rev. B 42, 7967 
(1990). 

12 Y. M. Vilk, Phys. Rev. B 55, 3870 (1997); Y. M. Vilk and 
A.-M. S. Tremblay, J. Phys. I France 7, 1309 (1997). 



1 M. Tinkham, Introduction to Superconductivity, 2nd ed. 
(McGraw-Hill, 1996). 

14 For similar considerations in the case of magnetism, see G. 
G. Lonzarich and L. Taillefer, J. Phys. C 18, 4339 (1985). 

15 We assume that the inclusion of pair breaking (represented 
by x) leads to a functional form for D equivalent to Eq. 4. 

16 The Kramers-Kronig is performed by linearizing ImE over 
small intervals, each of which is then analytically inte- 
grated and summed. 

17 W. J. Skocpol and M. Tinkham, Rep. Prog. Phys. 38, 1049 
(1975). 

18 N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 
J.-B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy 
and L. Taillefer, Nature 447, 565 (2007). 

19 P. A. Lee, Rep. Prog. Phys. 71, 012501 (2008). 

20 M. R. Norman, A. H. MacDonald and H. Akera, Phys. 
Rev. B 51, 5927 (1995). 

21 K.-T. Chen and P. A. Lee, Phys. Rev. B 79, 180510(R) 
(2009). 

22 A. J. Millis and M. R. Norman, Phys. Rev. B 76, 
220503(R) (2007). 

23 J. Orenstein, J. Corson, S. Oh and J. Eckstein, Ann. Phys. 



15, 596 (2006). 78, 17 (2006). 

P. A. Lee, N. Nagaosa and X.-G. Wen, Rev. Mod. Phys. 



